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Abstract 

We review and analyze the hybrid quantum-classical NMR computing methodology 
referred to as Type-II quantum computing. We show that all such algorithms consid- 
ered so far within this paradigm are equivalent to some classical lattice-Boltzmann 
scheme. We derive a sufficient and necessary constraint on the unitary operator rep- 
resenting the quantum mechanical part of the computation which ensures that the 
model reproduces the Boltzmann approximation of a lattice-gas model satisfying 
semi-detailed balance. Models which do not satisfy this constraint represent new 
lattice-Boltzmann schemes which cannot be formulated as the average over some 
underlying lattice gas. We close the paper with some discussion of the strengths, 
weaknesses and possible future direction of Type-II quantum computing. 
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1 Introduction 



The existence of quantum computers with thousands of error resistant entangled 
qubits is many years, if not many decades, in the future. However, the success of 
liquid state nuclear magnetic resonance spectroscopy (NMR) quantum computation 
means that few (~ 10) entangled qubit quantum computers are available now [1,2]. 
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This approach exploits the fact that the spin active nuclei in the molecules of a liquid 
are decoupled from each other and so each nuclear spin may act as a separate quantum 
bit, and each molecule as a separate few-qubit quantum computer. The states of these 
nuclear spins may be addressed by radio-frequency pulses. The existence of Avagadro's 
number of such spins enables the ensemble-averaged results of quantum computations 
to be obtained via thermodynamic averaging. 

The resonant frequency depends on the local magnetic field, and so placing the fluid 
in a non- uniform field enables one to separately address different regions of fluid [3] . 
Each region of fluid then acts as a separate ensemble of quantum computers, and 
the idea of using a large array of these computers was first proposed in [4]. Quantum 
coherence is maintained within each node on the array, and the arrays are measured 
after a few quantum coherent operations and then classical information representing 
the ensemble average of ~ 10^® measurements is exchanged between nodes, providing 
the initial state for the next coherent quantum computations on each node. This 
situation is illustrated in Figure 1. 

Type II quantum computers were introduced in the context of physical modelling, 
and all subsequent work has been focussed towards implementation of some lattice 
based algorithm for (usually classical) physics [5] . Algorithms have been proposed for 
fluid dynamics [6,7,4], for the diffusion equation [8], the Burgers equation [9], for the 
nonlinear Schroedinger and Korteweg-de Vries equations in one-dimensions [10], and 
for MHD turbulence in one dimension [11]. Some of these algorithms have also been 
implemented on NMR quantum computers [12,3,13]. 



2 Type II Quantum Lattice Gases 

The lattice structure of the array of quantum processors and the similarity of the 
two step quantum computation/classical communication are very similar to the col- 
lision/propagation evolution of lattice-gas and lattice-Boltzmann algorithms [14,15]. 
Because of this similarity Type-II algorithms are described in the literature as quan- 
tum lattice-gases. One must distinguish these Type-II models from the quantum lat- 
tice gases for the Dirac and Schroedinger equation, designed to realize exponential 
improvements in performance on fully coherent quantum computers [16,17]. Type 
II quantum lattice-gases perform the collision step using quantum computation, but 
then sample the results of those computations before performing the subsequent prop- 
agation. Type II simulations rely on ensemble averaging, and these simulations exploit 
the natural existence of an ensemble within the NMR quantum computer to estimate 
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Fig. 1. A single step of a Type-II quantum lattice-gas computation. Thick wires indicate 
classical registers. Thin wires indicate quantum registers. The computation proceeds by 
initializing a set of parallel quantum computations (the boxes "Q") representing a superpo- 
sition over the states of a given site of a classical lattice gas. The collision step is performed 
by a set of independent ensembles of quantum computations, one ensemble for each lattice 
site. The results of the measurements are a set of single particle occupation probabilities 
for the lattice vectors at each site, (classical information), are input to a classical computer 
(the box "C"), which performs the propagation step of the lattice gas algorithm. 

the average values of the qubits at a site. Note carefully, however, that only infor- 
mation about classical averages is communicated between sites. Information about 
two-particle and higher correlations is not preserved by this measurement and com- 
munication step, and so the results of a Type II quantum lattice gas simulation are 
equivalent to the results of a classical lattice Boltzmann simulation. Indeed, they 
most closely resemble the early lattice-Boltzmann algorithms which used ensemble 
averaged lattice gas collision rules at each site [18]. 

All of the models implemented so far on NMR quantum computers have two qubits 
per site in one dimension. The state of the model at a site is given by a pair of single 
particle distribution functions, (/i, /2), such that the probability that bit one is one 
is /i and the probability that bit two is one is /2. Each separate quantum processor 
therefore consists of a pair of coherent qubits. The state of these qubits is initialized 
using (/i,/2) as follows: 

m = 0)) = v/(i - /i)(i - /2) |oo)+v'(i - /i)/2 |oi)+v/(/i(i - /2) |io)+/Ai^|ii) 

(1) 
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The probabilities of the four possible lattice gas states at a site, {00,01, 10, 11} are 
n = ((1 - /i)(l - /s), (1 - /2)/i, /2(1 - /i), /1/2) in the Boltzmann approximation. 
The diagonal components of the density matrix \ip{t = 0)) {i/jit = 0)| are equal to 
these probabilities. The collision step of a stochastic lattice gas model is given by 
the multiplication of the vector n by a matrix. A, satisfying normalization and semi- 
detailed balance (a doubly stochastic or complete matrix). In our Type-II simulation 
we wish to reproduce the effect of this matrix multiplication by conjugating the 
density matrix by the unitary operator corresponding to the quantum mechanical 
part of our Type-II algorithm. The effect of this conjugation on the diagonal elements 
of the density matrix must be the same as the effect of matrix multiplication of vector 
n by the matrix A. 

Now consider a Type-II quantum computer with b qubits per node. At the beginning 
of each simulation step the qubits are initialized with a density matrix whose entries 
correspond to the probabilities of the qubit state in the molecular chaos approxi- 
mation. After each local unitary update, the qubits are measured and because the 
existence of an ensemble is presumed in all Type-II algorithms the probability that 
qubit i is one is associated with the single particle distribution function /j(x) at site x 
and vector i. During the collision step the unitary action U is applied to this density 
matrix and it evolves to: 

Pmr ' ^ ^ ^ UprnPrnr 
mr 

The density matrix formulation of quantum mechanics contains within it classical 
states, real linear convex combinations of classical states representing classical prob- 
ability distributions, pure quantum states and mixed quantum states. For this rea- 
son Type-II algorithms may be constructed which reproduce deterministic reversible 
lattice gas rules, stochastic lattice gas rules, and new models which have collision 
operators which do not correspond to the classical average over some lattice gas. In 
all cases, because the measurement and averaging step precedes the classical com- 
munication step, the overall algorithm is equivalent to a classical lattice-Boltzmann 
scheme with collision operator given by (2). 

Taking our update for the density matrix and restricting it to an update of the 
diagonal components of the density matrix: 

(3) 

Where there is no implicit sum on repeated indices. Separating the action on the 
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(2) 



Ppp ~ 51/ UprnPrnr U^ 
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diagonal and off diagonal components of the density matrix: 



Ppp 5Z ^pmUp^Pmm + UprnUp^Pmr (4) 



If the second term here is zero then the action of the unitary operator U on the 
diagonal components of p is identical to the action of a matrix A whose components 
are the modulus squared of the components of U. 

■^pm — UpmUp^ (5) 

The unitarity of U implies that the matrix A will obey normalization and semi- 
detailed balance. 

It remains to formulate a sufficient constraint that the second term above is zero. A 
sufficient condition for the Type-II algorithm to reproduce the action of a stochastic 
lattice-gas collision operator whose average is the matrix A is that: 

Re[Up^U;,]=0 ym<ryp (6) 



One realization of a Type-II algorithm on an NMR quantum computer implemented 
a simple stochastic lattice-gas algorithm for the diffusion equation in one- dimension, 
where the matrix A is: 



A=- 
2 



2 0^ 
110 
110 



(7) 
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In [3,12] this model is implemented using a unitary matrix U: 



U 



10 

Oif ii^O 
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It is clear by inspection that the modulus squared of the components of U give the 
components of A. The condition (6) is also satisfied. 
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If a Type-II algorithm does not obey constraint (6) the action of U on the diagonal 
components of the density matrix may not be described by a doubly stochastic matrix. 
The new values of the diagonal components will be modified by terms involving the off 
diagonal components of the density matrix, which in turn will involve the square roots 
of the single particle distribution functions. The action on the diagonal components 
of an arbitrary density matrix will not be linear and so (6) is a necessary as well as 
sufficient condition. Models violating (6) fall outside the class of lattice Boltzmann 
models which may be considered as averages over some underlying lattice-gas model 
in the Boltzmann approximation, and a detailed analysis of such lattice-Boltzmann 
models remains an open problem. An example of such a model is the lattice-gas model 
for the Burgers equation [19], which does not satisfy semi-detailed balance. A Type-II 
algorithm has been proposed for this model, and implemented in NMR [9,13]. 

In the case where the ensemble of molecular quantum computers at each node is 
initialized, not in the pure state (1), but in the completely mixed state with a density 
matrix with diagonal components equal to n and off diagonal components which are 
all zero, any unitary action reproduces the collision operator of a stochastic lattice gas 
satisfying detailed balance and normalization. In this case the Type-II computation is 
reduced to simple molecular computing of a Boltzmann approximation to a lattice-gas 
evolution, and no novel models are possible. 



3 Conclusions 

The characterization of current Type-II simulations as equivalent to classical lattice- 
Boltzmann methods illustrates the limited utility of Type-II quantum computation 
as presently defined. These algorithms provide at best a constant speedup and require 
an ensemble in order to accurately sample the classical information available through 
measurement which is exchanged between quantum processors. While applying tech- 
niques from medical imaging may make the realization of very large (> 2048^) grids 
possible, such techniques require a three dimensional field gradient to be useful. While 
the implementation of such field gradients is non-trivial, such large grid sizes would 
challenge the largest lattice Boltzmann simulations on classical supercomputers to 
date. 

One of the strengths of Type II quantum computation is that it takes advantage of 
the ensemble inherent in NMR implementations. Distribution functions are obtained 
by averaging over the ~ 10^^ molecules in each small region [20]. Of course, this 
necessity to obtain the average values by sampling becomes an Achilles heel once one 
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considers more scalable architectures (such as superconducting or solid state optically 
addressed implementations) where the necessity for sampling requires multiplication 
of the hardware. Although the spirit of Type II computing is that many small quantum 
computers are likely to be available sooner than a single large machine, the more 
scalable architectures are unlikely to reproduce the huge ensemble freely available 
with hquid state NMR. 

From the above analysis we can see that unitary matrices obeying (6) implement a 
doubly-stochastic Markov matrix action on the diagonal components of the density 
matrix. The size of this Markov matrix grows exponentially with the number of qubits. 
It is natural to ask whether one may use these Type-II techniques to solve Markov 
chain problems out of the reach of classical computers. For example, if we could 
address and control a pseudo-pure state of an ensemble of 100 qubit systems we 
could address Markov problems of size 2^^^, which are out of reach of any conceivable 
classical digital computer. However, the limit on the size of Markov process we may 
implement is not given by the number of qubits but by our ability to read out the 
diagonal components of the density matrix by sampling from the ensemble. We require 
at least one member of the ensemble per state of the system in order to do this, 
which implies that for 100 qubits one requires 10^ moles of material in our ensemble. 
Unfortunately this is not feasible. Current liquid state NMR realizations utilize 10^^ 
molecules, implying a maximum size of problem of 2^*^, implying that the maximum 
Markov problem size occurs for such an ensemble with 60 qubits. Such problems are 
classically intractable above at most 2^°, implying that there is a range of problems 
which can be implemented on small NMR ensemble quantum computers which can 
coherently address 30 — 60 spins which are classically intractable for conventional 
digital computers. 

This observation reinforces the point that Type-II quantum computation is really an 
interesting mixture of quantum computing with classical molecular computing. The 
approach benefits as much from the presence of a large number of realizations of the 
computation in a molecular ensemble as it does from the exponential growth of the 
dimension of the Hilbert space of the individual quantum processors. 

This material is based on work supported in part by the U.S. Army Research Labora- 
tory and the U.S. Army Research Office under grant number W911NF-04-1-0334, in 
part by the U.S. Air Force Office of Scientific Research under grant number FA9550- 
04-1-0176, and in part by the DARPA QuIST Program under AFOSR grant number 
F49620-01-1-0566. 
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